Nonhelical mean-field dynamos in a sheared turbulence 
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Mechanisms of nonhehcal large-scale dynamos (shear-current dynamo and effect of homogeneous 
kinetic helicity fluctuations with zero mean) in a homogeneous turbulence with large-scale shear 
are discussed. We have found that the shear-current dynamo can act even in random flows with 
small Reynolds numbers. However, in this case mean-field dynamo requires small magnetic Prandtl 
numbers (i.e., Pm < Pm'^'^ < f). The threshold in the magnetic Prandtl number, Pm'^'^ = 0.24, is 
determined using second order correlation approximation (or first-order smoothing approximation) 
for a background random flow with a scale-dependent viscous correlation time Tc — {uk^)~^ (where 
1/ is the kinematic viscosity of the fluid and k is the wave number) . For turbulent flows with large 
Reynolds numbers shear-current dynamo occurs for arbitrary magnetic Prandtl numbers. This 
dynamo effect represents a very generic mechanism for generating large-scale magnetic fields in a 
broad class of astrophysical turbulent systems with large-scale shear. On the other hand, mean- 
field dynamo due to homogeneous kinetic helicity fluctuations alone in a sheared turbulence is not 
realistic for a broad class of astrophysical systems because it requires a very specific random forcing 
of kinetic helicity fluctuations that contains, e.g., low- frequency oscillations. 
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I. INTRODUCTION 

An origin of solar, stellar and galactic large-scale mag- 
netic fields is related to a mean-field dynamo. This dy- 
namo can be driven by the joint action of small-scale 
turbulent flows with a nonzero mean kinetic helicity and 
large-scale differential rotation (see, e.g., Moffatt 1978; 
Parker 1979; Krause and Radler 1980; Zeldovich et al. 
1983; Ruzmaikin et al. 1988; Ossendrijver 2003; Riidiger 
and HoUerbach 2004; Brandenburg and Subramanian 
2005). 

Recent numerical experiments by Yousef et al. (2008) 
have demonstrated existence of a nonhelical large-scale 
dynamo in a turbulence with superimposed linear shear 
in elongated shearing boxes whereby mean a effect van- 
ishes. The exponential growth of magnetic field has been 
found at scales which are much larger than the outer 
scale of the turbulence. An earlier indications of non- 
helical turbulence amplifying large-scale magnetic field in 
the presence of a large-scale shear associated with mean 
fiows has been found by Brandenburg (2005) and Bran- 
denburg et al. (2005) in numerical experiments that used 
constant-in-time sinusoidal forcing functions. This im- 
plies that the amplification effect in a sheared nonhelical 
turbulence appears to be numerically robust. Note also 
that numerical experiments with Taylor-Green forcing is 
another example of a mean-field dynamo produced by a 
combined effect of a nonhelical turbulence and a compli- 
cated large-scale flow (Ponty et al. 2005). 

One of the possible mechanism of the nonhelical large- 
scale dynamo in a homogeneous sheared turbulence is a 
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shear-current dynamo effect (see Rogachevskii and Klee- 
orin 2003; 2004; 2007; Rogachevskii et al. 2006a, 2006b). 
The physics of this phenomenon is following. Upward and 
downward turbulent eddies result in deformations of the 
original nonuniform magnetic field lines. In a turbulence 
with a large-scale shear the inhomogeneity of the origi- 
nal mean magnetic field breaks a symmetry between the 
influence of the upward and downward turbulent eddies 
on the mean magnetic field. This causes the mean elec- 
tric current along the mean magnetic field and results in 
systematic amplification of the large-scale magnetic field. 

The shear-current dynamo has been previously stud- 
ied for large Reynolds numbers using the spectral tau- 
approximation (see Rogachevskii and Kleeorin 2003; 
2004; 2007). On the other hand, in a random flow 
with small Reynolds numbers the dynamo action in 
nonhelical shear flows has not yet been found in the 
framework of the second order correlation approxima- 
tion or flrst-order smoothing approximation (see Radler 
and Stepanov 2006; Riidiger and Kichatinov 2006). How 
generic is the latter result? In particular, how the dy- 
namo action may depend on the spectral properties of a 
random flow with small Reynolds numbers? One of the 
goals of this study is to revise this problem for the case 
of a random flow with small Reynolds numbers. 

Another effect that might explain the large-scale dy- 
namo in a sheared turbulence with a zero mean kinetic 
helicity is associated with kinetic helicity fluctuations. 
Dynamics of large-scale magnetic field in the presence of 
kinetic helicity fluctuations with a zero mean in a shear- 
free turbulence has been studied for the first time by 
Kraichnan (1976). This problem is formulated in the 
following way. Let us consider a small-scale turbulence 
produced by a random forcing F^"' located in small scales 
li> < ^turb < ^0 (and Tj, < Tturb < To), whilc in larger 
scales lo ^ I -^ Ix (and ^o *^ ''' ^ ^x) there are kinetic 
helicity fluctuations (or a fluctuations) with a zero mean 



produced by a random forcing F^''^'. 



II. THE SHEAR-CURRENT EFFECT 



The mean-field effects occur at very large scales L ^ l^ 
(and times r^ ^ r^), where the mean kinetic helicity 
is zero. All mean quantities are determined by double 
averaging over velocity fluctuations, (...), and over ki- 
netic helicity fluctuations (•■•)^"^ [see detailed discussion 
by Sokolov (1997) about various mathematical aspects 
of this problem] . Numerical simulations of the magnetic 
field evolution in accretion discs by Vishniac and Bran- 
denburg (1997) have demonstrated that kinetic helicity 
fluctuations with a zero mean can result in generation 
of large-scale magnetic field (see also Brandenburg et al. 
2008). 

Let us discuss theoretical aspects of this problem. In a 
shear-free turbulence kinetic helicity fluctuations cause 
two effects: (i) a negative contribution to the turbu- 
lent magnetic diffusion, 77^") — — t^((5^)("^; and (ii) 
a large-scale drift velocity of the mean magnetic field, 
V(") ex r^ V(a2)(") (see Kraichnan 1976; Moffatt 1978). 
In a turbulence with large-scale shear, inhomogeneous 
kinetic helicity fluctuations can produce a mean-field dy- 
namo (Silant'ev 2000). Indeed, a combined effect of the 
inhomogeneous fluctuations and large-scale shear super- 
imposed on turbulence, produces a nonzero mean al- 
pha effect: a(^'") oc -t^SV ,{a'^)'^'^\ while (a)'") = 0. 
Here the mean vorticity due to the large-scale shear is 
W = Sbz. The equation for a^-^'") has been derived us- 
ing the second order correlation approximation and the 
spectral tau-approximation (see for details, Kleeorin and 
Rogachevskii 2008). The large-scale shear and the mean 
alpha effect can result in the mean-field dynamo that acts 
similarly to the aQ-dynamo. 

Using phenomenological arguments, Proctor (2007) 
has suggested that homogeneous kinetic helicity fluctu- 
ations in a homogeneous turbulence with a large-scale 
shear may generate a large-scale magnetic field. Such 
possibility for a large-scale dynamo has been recently ex- 
amined by Kleeorin and Rogachevskii (2008) using the 
second order correlation approximation and the spectral 
tau-approximation. This study has not found large-scale 
dynamo produced by homogeneous kinetic helicity fluc- 
tuations alone with zero mean value in a sheared homo- 
geneous turbulence. However, how generic is the latter 
statement? One of the goals of this study is to revise 
this problem. We have demonstrated that only for a spe- 
cific random forcing of kinetic helicity fiuctuations that 
also contains low-frequency oscillations, the large-scale 
dynamo in a homogeneous turbulence with a large-scale 
shear might be possible. 

This paper is organized as follows. In Sec. II we investi- 
gate shear-current dynamo for a random flow with small 
Reynolds numbers and different spectral properties. In 
Sec. HI we study the effect of homogeneous kinetic helic- 
ity fluctuations with a zero mean in a sheared turbulence. 
In Sec. IV we draw concluding remarks. 



In order to study the shear-current effect in a random 
flow with small Reynolds numbers we use a second or- 
der correlation approximation (SOCA). This approxima- 
tion is valid only for small hydrodynamic Reynolds num- 
bers. Even in a highly conductivity limit (large magnetic 
Reynolds numbers), SOCA can be valid only for small 
Strouhal numbers (i.e., for very short correlation time). 

We use equation of motion and induction equation for 
fluctuations of velocity and magnetic fields, exclude the 
pressure term from the equation of motion by calcula- 
tion Vx(Vxu), where u are velocity fluctuations. We 
rewrite the obtained equation and the induction equation 
in a Fourier space and apply the two-scale approach (i.e., 
we use large-scale and small-scale variables). We neglect 
nonlinear terms but keep molecular dissipative terms in 
the equations for fluctuations of velocity and magnetic 
fields. We seek for a solution of the obtained equations 
for fluctuations of velocity, u, and magnetic, b, fields as 
an expansion for a weak velocity shear: 



b(k,cj) = b(°)(k,cj)+b(l)(k,w) + 
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Here U and B are the mean velocity and magnetic 
fields, k and lo are the wave number and frequency, 
Gi,{k, oj) = (z^fc^ — iuj)^^ and G^(fc, uj) — [rik"^ — iuj)~^, rj 
is the magnetic diffusion coefficient due to electrical con- 
ductivity of the fluid, and v is the kinematic viscosity of 
the fluid. For derivation of Eqs. I©-© we use an identity 

j C7,(Q) 6„(k - Q) rfQ = z(Vpi7,) ^ , 

that is valid at least for a linear velocity field. Equa- 
tions (121)-® coincide with that derived by Radler and 
Stepanov (2006). These equations allow us to determine 



the cross- helicity tensor gmn — {um bn') + {um bn ) and 
the contributions, £^ = Simn J ffmn(k, a->) dkdo;, to the 
electromotive force caused by sheared turbulence. For 
the integration in k-space we have to specify a model 
for the background shear- free turbulence (with B = 0), 
which is determined by equation: 



[Ui Uj 



.(0)(k,C^) = (u2)(") 



P,,{k)E{k) 



87r2 fc2 r^ (^2 



') 
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where E{k) is the energy spectrum (e.g., a power-law 
spectrum), Tc is the correlation time, Pij{k) — 5ij — 
kikj/k"^ and Sij is the Kronecker tensor. This model cor- 
responds to the correlation function: {ui{t)uj{t + r)) ex 
exp(— t/tc). Straightforward calculations yields the con- 
tributions to the electromotive force caused by sheared 
turbulence: 



AS) 



£l°' = lt,[A^e,pk{dU)p^{dB)qk+A2Wk{dB),k 
+A3Jfc(5C7),fe+A4(WxJ),], 



(7) 



where {dU),j = {V,Uj + VjC/j)/2, W = VxU is the 
mean vorticity, Iq is the maximum scale of turbulent mo- 
tions (the energy containing scale), J = VxB is the 
mean electric current, and the coefRcients An are given 
in Appendix. The equation for the evolution of the mean 
magnetic field, B = [Bx{z),By{z),Q), reads 



dt 

dB^ 

dt 



= -a,SllB'^ + iii + 7^,)B'^, (8) 

= SBx + {7^ + r,^)B'' (9) 



where we use linear velocity shear U = (0, Sx, 0), B'/ = 
d^Bi/dz^, rj^ c>c tq (u^)^"-' is the turbulent magnetic dif- 
fusion coefficient, tq = lo/\/{u^)^^\ and we neglect small 
contributions to the coefficient of turbulent magnetic dif- 
fusion caused by the shear motions because we consider 
a small shear, Stq <C 1. The coefficient a^ entering in 
Eq. ^ is given by 

^ f[h- uk^h + Vk^h - h)] E{k) fc2 dk , 

(10) 
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the functions In{k) for Tc = l/{vk'^) are given in Ap- 
pendix. Using the explicit form of the functions /„(fc), 
we obtain the following expression for the coefficient a^ : 



60 (to v)- 



■ Pm (1 - 4Pm - Pm^ 



E{k) 



dk. 
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where Pm = v/rj is the magnetic Prandtl number. The 
solution of Eqs. ([8]) and ^ we seek for in the form ex 
exp(7^ t + iKz z), where the growth rate, 7^ , of the mean 
magnetic field is given by 



and cr^ > when Pm < 0.24. 

In the present study we use the SOCA procedure that 
is valid only for Re <C 1. It follows from Eqs. (fTT|) 
and p^ that for Re <C 1, the dynamo instability due to 
the shear-current effect occurs when Pm < 0.24 (i.e., for 
small magnetic Prandtl numbers). This result has been 
obtained for a model of the background shear-free tur- 
bulence determined by Eq. © with Tc = l/{vk'^). Note 
that Radler and Stepanov (2006) used a model of the 
background shear-free turbulence with a constant scale- 
independent correlation time Tc- A possibility for the 
shear-current dynamo for small magnetic Prandtl num- 
bers in the case of Re <C 1 has been pointed out by 
Riidiger (2007), although this was not explicitly men- 
tioned in his previous study using the SOCA procedure 
and a more simple model for the background shear-free 
turbulence: {uiU.j)^°H\i,uj) ex {u^)^°'i P,j{k) E{k) 5{uj) 
(see Riidiger and Kitchatinov 2006). For turbulent flows 
with large Reynolds numbers shear-current dynamo oc- 
curs for arbitrary magnetic Prandtl numbers (see Ro- 
gachevskn and Kleeorin 2003; 2004; 2007). 



III. EFFECT OF KINETIC HELICITY 
FLUCTUATIONS 

In order to study effect of kinetic helicity fluctuations 
with a zero mean on large-scale dynamo we use a second 
order correlation approximation. This procedure yields 
the equation for the evolution of the magnetic field B: 

— = Vx(^5B + VxB -(77 + 77^ )Jj+B^, (13) 

where J = VxB is the electric current, B^ are the non- 
linear terms, V -I- u is the total velocity and (u) — 0. In 
this section we do not consider the shear-current effect. 

In the scales Iq -^ I -^ l^ there are fluctuations of 
a. Let us consider homogeneous kinetic helicity fluctu- 
ations. In order to derive equation for the the mean 
magnetic field B = (B)("\ we determine the contri- 
bution to the mean electromotive force caused by the 
sheared turbulence and the kinetic helicity fluctuations. 



(S.Q) 



£}.^^°'> = (aBj)^°'\ To this end we use Eq. ^ in which 
we neglect the nonlinear terms B^. Solving this equa- 
tion in a Fourier space we determine the magnetic field 
-Bj(k,a;), where the wave vector k and the frequency lo 
are located in the spatial scales Iq -^ I -^ l^ and in the 
time scales tq -^ t ^ t^. Multiplying the magnetic field 
ByCkjLu) by a and averaging over kinetic helicity fluctu- 
ations we determine £■, 
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where GrikjUj) = [(77 -I- r]^)k'^ — Jw]"""", fa{k,uj) — 
(a(w)a(— w))'^"\ U = {0,Sx,0) is the background shear 
velocity, and J = VxB is the mean electric current. 
We assume that the mean magnetic field has the form: 



B = {Bx{z), By{z),0), and neglect small contributions 

~ 0{to/t-^) to the mean electromotive force £y '" . 

We use the following model for the spectral function 
fa{k,uj): 



faik^Uj) 



(a^\{») 



E^{k) 



TTT^iiO^+T^ ) 



(15) 



This model corresponds to the following correlation func- 
tion {a{t)a{t + t))("^ cx exp(— t/t^). In earher stud- 
ies by Kleeorin and Rogachevskii (2008), a more simple 
model for the spectral function has been used: /a(fc, ^) = 
(a2)(")^„(fc)5(tj). 

The contribution to the mean electromotive force 
caused by the sheared turbulence and the kinetic helic- 
ity fluctuations is given by fj '"' = -CTq S (q;^)^") t^ By, 
where the parameter 



E^{k) 



[l + rx(?7 + r?^)fc2] 



■dk>0 



(16) 



Here we use an identity / G^ Gq G* duj = tt/ [a [rj k^ + 
a)^], where Ga{k,uj) = (a — ico)^^ with o = r^^. The 
equation for the evolution of the mean magnetic field, 
B== {B^{z),By{z),0), reads 



dt 

dt 



ao^S{a'Y"K^B'' + f^,B':, (17) 



SB^ + fj^B'^ 



(18) 



where fjj, = 77 + 77^ + rjl^' . Here we neglect small contri- 
butions to the coefficient of turbulent magnetic diffusion 
caused by the shear motions because Stq ^ 1. Note that 
for enough general model pS]) of the spectral function 
fa{uj,k), the parameter CTq, is always positive. It is also 
positive when r^ = (i^fc^)"^ (see the spectral model used 
in Sect. 2). This implies that homogeneous kinetic he- 
licity fluctuations alone with zero mean value for general 
model (jlSp in a sheared homogeneous turbulence cannot 
cause a large-scale dynamo [see Eqs. (|T7| and ([18])]. 

However, for a specific random forcing of kinetic helic- 
ity fluctuations that also contains low-frequency oscilla- 
tions, e.g., 

{a{t)a{t + t))(") cx exp(-T/rx) cos(ww'r) , (19) 

there is a possibility for a large-scale dynamo action due 
to homogeneous kinetic helicity fluctuations in a sheared 
homogeneous turbulence. In this case the spectral func- 
tion fa{k,uj) is given by: 

-2\(a) Ea{k) 



/„(fc,c.) = (a2)W 



27r 



Lj2 + [a*Y 



with a = T^ 



and the parameter a a, 



{[l + Tx(r7 + r7^)fc2]2 + (^^^^)2}; 



■Eaik)dk 



(20) 



is negative when uj^jT^ > 1 + t^ {ij + r]^) k"^ . This implies 
that for model p^ a large-scale dynamo due to homoge- 
neous kinetic helicity fluctuations in a sheared homoge- 
neous turbulence can occurs. However, this model for the 
function /^(fc, to) of kinetic helicity fluctuations seems to 
be not realistic. 



IV. CONCLUSIONS 

Two types of nonhelical large-scale dynamos due to 
shear-current effect and homogeneous kinetic helicity 
ffuctuations with zero mean in a sheared turbulence are 
investigated using a second order correlation approxima- 
tion. The mechanism for the shear-current dynamo is fol- 
lowing. The large-scale velocity shear creates anisotropy 
of turbulence that produces a contribution to the electro- 
motive force, WxJ, caused by the shear. Joint effects 
of the electromotive force WxJ and stretching of the 
mean magnetic field due to the large-scale shear motions 
cause the shear-current dynamo instability. This effect 
occurs even for small Reynolds numbers. However, the 
dynamo instability in this case requires small magnetic 
Prandtl numbers (Pm < 0.24). This dynamo threshold is 
found for a model of a random ffow with the correlation 
time Tc — {vk"^)^^ . The shear-current dynamo for large 
Reynolds numbers is independent of magnetic Prandtl 
numbers. 

Another possible mechanism for the nonhelical large- 
scale dynamo is associated with homogeneous kinetic he- 
licity fluctuations in a sheared turbulence. However, this 
kind of mean-field dynamo is not universal and can occur 
only for a specific random forcing of kinetic helicity fiuc- 
tuations that contains, e.g., low- frequency oscillations. 

The discussed effects in this study might be impor- 
tant in a broad class of astrophysical flows. For instance, 
sheared turbulence is a universal feature in astrophysi- 
cal flows, e.g., in stellar interiors, accretion disks, irreg- 
ular galaxies (Balbus and Hawley 1998; Chyzy et al. 
2000; Ossendrijver 2003; Brandenburg and Subramanian 
2005; Donati et al. 2005; Gaensler et al. 2005), and in 
liquid- metal laboratory dynamo experiments (see, e.g., 
Monchaux et al. 2007). 

Non-symmetrical explosions of supernova may produce 
fluctuations of kinetic helicity located in larger scales 
than small-scale turbulence existing in convective zones 
inside stars. On the other hand, the shear-current dy- 
namo acts together with the a-shear dynamo. The shear- 
current effect does not quenched (see Rogachevskii and 
Kleeorin 2004; Rogachevskii et al. 2006b) contrary to the 
quenching of the nonlinear a effect, the turbulent mag- 
netic diffusion, etc. This implies that the shear-current 
dynamo might be the only surviving effect, which can ex- 
plain the origin of large-scale magnetic fields in sheared 
astrophysical turbulence. 
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APPENDIX A: COEFFICIENTS A„ AND 
FUNCTIONS /„(fc) 

The coefficients A„ entering in Eq. ([7]) are given by 
Ai = -^ [[2h-2h + 5Is + 2iyk\2l3+l6) 

OUTT J 



A. 



-^2r]k^{2l2 ~ I7)] E{k) k^ dk , 

-^ [[6l4 + 5Is-6j^k^h + 2rjk^{3Ii-3l2 
bvn J 



-2/7)] E{k) k^ dk , 

^3 = 7^ f[5h-2h-2iyk\h+Ie) 
bOn J 



-2rik^{Ii 
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24 TT 
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where the functions In{k) for t^^ — vk^ are given by 
h{k) = yC^G^G: da; =2^2(^ + ^)2 ^8 ' 
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